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C^ , Abstract. We present integral field spectroscopy in the near infrared (NIR) of He 2-10 and NGC 5253, two well 

known nearby dwarf irregular galaxies showing high star-formation rates. Our data provide an unprecedented 
*ivj , detailed view of the interstellar medium and star formation in these galaxies, allowing us to obtain spatially 

^ . resolved information from the NIR emission and absorption line tracers. We study the spatial distribution and 

\j^ ' kinematics of different components of the interstellar medium (ISM) mostly through the Bracket series lines, the 

V^ , molecular hydrogen spectrum, [Fell] emission, and CO absorptions. Although the ISM is mostly photo-excited, 

V^ ' as derived by the [FeII]/Br7 and H2 line ratios, some regions corresponding to non-thermal radio sources show a 

CO ' [FeII]/Br7 excess due to a significant contribution of SN driven shocks. In He 2-10 we find that the molecular gas 

r^^ . clouds, as traced by CO(2-l) and H2 infrared line, show consistent morphologies and velocities when studied with 

f^ ' the two different tracers. Moreover, there is a clear association with the youngest super star clusters as traced by 

f^ I the ionized gas. In the same galaxy we observe a cavity depleted of gas, which is surrounded by some of the most 

active regions of star formation, that we interpret as a signature of feedback-induced star formation from older 

episodes of star formation. Finally, we measured high turbulence in the ISM of both galaxies, a ~ 30 — 80 km/s, 

driven by the high star-formation activity. 



?H ' Key words. Galaxies: dwarf - Galaxies: Individual: He 2-10 Galaxies: Individual: NGC 5253 - Galaxies: ISM 



Galaxies: starburst 

1. Introduction evolution. Starburst galaxies are complex systems, where 

gas is converted into stars on a relatively short timescale 
Detailed observations of intense star formation and mass ^-^^ ^^^^^ efHciency: the theoretical progress in this area 
assembly are a key ingredient in probing the role of various ^^ ^^-^^ ^^^^^^ dependent on highly simplified recipes for the 
mechanisms m shapmg the galaxies and moderating their physical mechanisms that drive the formation of new stars 
■ ^ and their interplay with the interstellar medium (ISM). 

■"^ ^ ' ■ Building a complete theoretical picture is important to 
Based on observations obtained at the ESO-VLT under 



programs 074. A- 9011 and 075.B-0648 



fully understand the mechanism responsible for the for- 
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mation and evolution of star-forming galaxies at high red- 
shift (z > 1), where the bulk of star formation in the 
Universe took place (e.g. Rudnick et al. I2006p . Important 
progress has been made in recent years by direct observa- 
tions of distant galaxies, but at those distances even the 
largest sources are just a few arcsec in size (e.g. Cresci et 
al. I2006p . preventing us from distinguishing the different 
mechanisms at work. Although perfect analogs of high- 
z galaxies are not available locally owing to the different 
conditions in the ISM, gas content, accretion rate, and 
chemical properties, studying star-forming galaxies at a 
much higher physical resolution is critical to investigate 
some key properties that are virtually inaccessible at high 
redshift. Nearby starburst systems have already helped in 
the past decades to understand some of the processes that 
regulate and trigger star formation on galactic scales, e.g. 
the interaction with close companions or large molecular 
clouds, and internal mechanisms such as feedback from 
previous generations of massive stars. Blue compact dwarf 
galaxies are particularly suited for this task, as they show 
typically sub-solar abundances, thus providing the oppor- 
tunity to observe star formation in a relatively pristine 
environment and with less contamination by the underly- 
ing evolved stellar population than in larger galaxies. 

The population of young super massive stellar clusters 
(SSCs) that accounts for a large fraction of the star for- 
mation of the entire galaxy was identified in starburst sys- 
tems both locally and at high-z (e.g. de Grijs et al. :20031 
Elmegreen et al. l2005p . These clusters could be the progen- 
itors of the old globular clusters observed in most galax- 
ies today (de Grijs & Parmentier 12007^ . thus providing 
further evidence for the link between the starburst phe- 
nomenon and the formation of galaxies. Moreover, star- 
forming regions with properties consistent with the more 
massive HII regions detected locally are observed up to 
z ~ 5 (Swinbank et al. I2009L and local dwarf galaxies 
have been identified as scaled down analogs of clumpy 
star- forming galaxies at high redshift (Elmegreen et al. 
I2009P . Dwarf galaxies indeed show high gas fraction, high 
turbulence, a large Jeans length compared to their size and 
a clumpy appearance, similar to what is observed in high- 
z forming disks (see e.g. Genzel et al. 12006120081 Cresci et 
al. 12009. Forster-Schreiber et al.HOMl)- In this picture, the 
clumpy morphology comes from gravitational instabilities 
in gas with high turbulent speed compared to the rota- 
tional one, and high gas mass fraction compared to stars. 
Although the mechanism providing the gas in the local 
and high-z galaxies might be different (mergers and inter- 
actions are required locally, while cosmological accretion 
through cold flows might be dominant at high redshift, see 
e.g. Dekel et al. l2009p . the physical process regulating the 
star formation in clusters are probably the same in the 
two environments. 

In the effort of understanding the physical processes at 
play in galactic-scale star formation, we thus undertook a 
survey campaign of nearby dwarf starburst galaxies us- 
ing integral field spectroscopy in the near infrared (NIR). 
This technique provides an unprecedented detailed view 



of the interstellar medium and star formation in these 
galaxies, allowing us to obtain spatially resolved informa- 
tion from the emission and absorption line tracers in the 
NIR. The results derived from the study of the nearby 
starburst galaxy II Zw 40 were presented in Vanzi et al. 
(|2008p . while in this paper we discuss the observations ob- 
tained on He 2-10 and NGC 5253, two very well known 
examples of nearby starbursts. Both galaxies are dwarf- 
irregular, with estimated distances of 9 Mpc for He 2-10 
(Vacca & Conti [T992ll and 3.3 Mpc for NGC 5253 (Gibson 
et al. I2000p : due to their relatively close distances and re- 
markable properties they have been studied in detail in 
the past. 

1.1. He 2-10 

He 2-10 has nearly solar metahcity {l2 + log{0 / H) = 8.93; 
Kobulnicky et al. I1999ap and hosts an intense episode of 
star formation as revealed by classical indicators such as 
the strong optical and infrared emission lines (Vacca & 
Conti ITM^ Vanzi & Rieke [TM7|) . the Wolf-Rayet fea- 
tures (Vacca & Conti I1992p . intense mid-IR (Sauvage 
et al. I1997P and far-IR continuum (Johansson I1987p . 
Some compact radio sources were observed by Johnson 
& Kobulnicky (|2003p . who characterized their spectra as 
thermal, measured their sizes as 2 — 4 pc and their masses 
in excess of 10^ Mq. A number of young super massive 
clusters in the star-forming region of the galaxy, with typ- 
ical masses up to 10^ Mq and ages of few Myr have been 
found in optical and UV studies (e.g. Johnson et al. 120001 
Vacca & Conti IHnil). However, these observations reveal 
clusters that are relatively unobscured and are proba- 
bly not the youngest star-forming regions. Cabanac et al. 
P005l hereafter C05) observed in the IR, to penetrate the 
obscuring material, and detected several compact sources 
that could be associated with the radio sources previously 
mentioned, finding that most of them are ultra-compact 
HII regions and two are possible SN remnants (see Fig. 

ED. 

Kobulnicky et al. (|1995p derived a total molecular gas 
mass of 1.6 • 10^ M© and an atomic gas mass of 1.9 • 
10® Mq through single-dish observations of CO(l-O) and 
C0(2-l) and VLA HI 21cm. Thanks to the interferometric 
resolution, they detected a CO plume extended about 30" 
to the southeast, and a similar but more extended feature 
in HI. They interpreted these features as a sign of a merger 
with a companion galaxy. This interpretation is supported 
by the detection of an extended component of CO (3-2) 
molecular gas to the northeast of the main body of the 
galaxy by Vanzi et al. (|2009p . interpreted as additional 
evidence of a detached cloud expelled in the merger that 
He 2-10 is experiencing. 

1.2. NGC 5253 

NGC 5253 is an irregular dwarf galaxy in the Centaurus 
Group and another interesting target to study young star- 
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bursts. It is indeed one of the closest starburst galaxies, 
with a heliocentric distance of only 3.3 Mpc (Gibson et 
al. I2000p . so that 1" corresponds to 16 parsec, and one of 
the youngest starburst galaxies known (Rieke et al. ll988p . 
as implied by the detection of spectral features arising 
from Wolf-Rayet stars (Schaerer et al. I1997P and from 
its almost entirely thermal radio spectrum with very lit- 
tle synchrotron emission from supernova remnants (Beck 
et al. I1996p . Its metalicity is sub-solar and about Zq/Q 
(Kobulnicky et al. 1999). 

At centimetric wavelengths, Turner et al. (|1998p and 
Turner et al. (|2000p have observed several nebulae, which 
they derive to be ionized by 200 — 1000 O stars, suggesting 
that very large clusters are the preferred mode of star 
formation in the central regions of the galaxy. Calzetti 
et al. (|1997p found that they span a range of masses 
and ages, and identified six main clusters, one of which 
powers a giant HII region that is the most prominent 
object in the galaxy (Vanzi & Sauvaee l2004p . Cresci et 
al. (|2005p detected more than 300 clusters combining 
high-resolution NIR and optical imaging, with derived 
ages spanning a range between 2 and 50 Myr and masses 
of up to 10*^ Mq. 

The paper is organized as follows: Sect. [5] provides the 
details of the observations, in Sect.[3]we describe in detail 
the results, in Sect. U we discuss a possible origin of the 
cavity detected in the ISM of He 2-10, while in Sect. [S]we 
summarize our results. 

2. Observations and data reduction 

The observations were obtained with the integral field 
spectrograph SINFONI at the ESO-VLT (Eisenhauer et 
al lMIHl Bonnet et al. IMH)) as part of the ESQ program 
075.B-0648 (He2-10 and NGC5253) and 074.A-9011 (Hc2- 
10). 

We used the 0.250" x 0.125" pixel scale, which provides 
a field of view of 8" x 8". Due to the large apparent sizes 
of these two galaxies, we were able to cover only part of 
the optical extent of the targets with SINFONI. For He 2- 
10 our pointing is centered on the bright central nucleus, 
generally referred to as region A in the literature, which is 
an arc of UV-bright super star clusters that dominate the 
emission of the galaxy (Vacca & Conti lTW^ . In NGC 5253 
our observations were centered at the brightest HII region, 
i.e. cluster 5 according to the notation of Calzetti et al. 
(fT997)l . 

For each galaxy we obtained two data cubes using 
the H and K gratings, providing one spectral resolution 
R=3000 and 4000 respectively. He2-10 was observed on 
2005 March 24 and 2005 April 3 for a total integration 
time on-source of 20' in both H and K band. NGC5253 
was observed on 2005 April 3, with an integration time of 
10' in both bands. 

The data were reduced with the SINFONI custom re- 
duction package SPRED (Abutcr et al. I2006P , which in- 
cludes all of the typical reduction steps applied to near-IR 



spectra with the additional routines necessary to recon- 
struct the data cube. After background subtraction, which 
was performed with off sky frames obtained interspersed 
with the on-source exposures, the data were flat-fielded 
and corrected for dead/hot pixels. Telluric correction and 
flux calibration were carried out using A- and B-type stars, 
after removing the stellar features. Residuals from the OH 
line emission were minimized with the methods outlined 
in Davies ([MITp . 

A crude estimate of the seeing limited spatial resolu- 
tion achieved was obtained with the standard star obser- 
vations. These show a FWHM ~ 0.6" for both galaxies, 
corresponding to a linear resolution of 26.2 and 9.6 pc at 
the distances assumed for He 2-10 and NGC 5253 respec- 
tively. 

An example of the obtained very rich spectra for the 
two galaxies is shown in Fig. [5] Our analysis is focused 
on the main visible features, mainly the emission lines of 
atomic hydrogen, the emission lines of molecular hydro- 
gen, the forbidden emission lines of [Fell], and the stellar 
absorption features of CO. We extracted spectral images 
centered on the main features detected, fitting each emis- 
sion or absorption line at each spaxel with an unresolved 
line profile (e.g. a sky line) convolved with a Gaussian, as 
well as a linear function to the line-free continuum. This 
fitting process gives the intrinsic velocity dispersion (i.e., 
instrumental broadening has been accounted for), the ab- 
solute velocity and the total fiux in the line at that loca- 
tion, and these values are used in the subsequent analy- 
sis, figures, and tables. The uncertainty of the fit is esti- 
mated with Monte Carlo techniques, by refitting the best- 
fit Gaussian with added noise of the same statistics as the 
data. This is done 100 times, and the standard deviations 
of the best-fit parameters are used as the uncertainties 
(see Cresci et al. IMM Forster Schreiber et al. 120051) . 

3. Analysis 

In Figs. [3] and 2] we show the line maps, extracted as ex- 
plained in Sect.m in Br7 2.167 /^m, [Fell] 1.643 ^m, H2 
2.121 /xm rest wavelength for both galaxies. The contours 
of the K-band continuum from the SINFONI data-cubes 
are drawn for comparison. Below we will mainly use these 
tracers to study the properties and dynamics of the dif- 
ferent phases of the ISM in the two galaxies: ionized gas 
through recombination lines of H and other species like 
Fe, warm molecular gas through H2, and dust extinction 
using ratios of well-modeled lines. 

3.1. Ionized hydrogen 

The ionized gas distribution, which we can best trace in 
our data through the Br7 emission, allows us to identify 
the main regions of star formation in the two galaxies. We 
refer to the first panel of Figs. [3] and S] for He 2-10 and 
NGC 5253 respectively. 

In He 2-10 three main regions are prominent in Br7. 
The central region, where the emission peaked, is also 
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Fig. 1. K-band continuum images of He 2-10 (left panel) and NGC 5253 (right panel), extracted from the SINFONI 
datacube. The Br7 contours are over-plotted as reference. The main clusters and star-forming regions in the field of 
view are identified with the naming convention from literature: for He 2-10, the clusters (Ll-6) from Cabanac et al. 
(|2005p and the radio knots (knot 1-5) from Kobulnicky & Johnson (|1999p are marked, while for NGC 5253 the clusters 
from Calzetti et al. (|1997p and Vanzi & Sauvage (|2004p . 
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Fig. 2. Near-infrared H and K spectrum extracted with a circular aperture of 0.375" radius around the Br7 peak of 
the two galaxies, with the main detected emission lines explicitly indicated, execept for the large number of lines from 
the Brackett series detected in the H-band. 



bright in Ha, which coincides with the region called L4 
by C05. This area is bright and extended in the radio 
continuum and is identified as knot 4 by Kobulnicky & 
Johnson p9991 hereafter KJ99). A second bright Br7 re- 
gion is seen to the southeast of the previous one. It is not 
bright in Ha, but it is detected in the I, K, and N band 
continuum (C05), probably due to higher extinction (see 
Sect. 13. 4|) . It coincides with source L5 of C05 and with the 
radio knot 5 of KJ99. Finally, a region to the northwest of 
the central one is also detected in Ha as a diffuse nebula, 
bright in L and 10 fim (knot LI in C05) and coincident 
with the radio sources knot 1 and 2 of K J99 (see Fig. [T]). 
The K continuum detected in SINFONI and shown 
in Figs. [1] and [3] in contours is identical to the K image 
obtained by CG5. The comparison of the Br7 morphology 
and the continuum shows a clear difference, with the max- 



imum of the K continuum located 0.7" to the NW of the 
Br7 one, and a plume of three clusters bright in contin- 
uum but undetected in Br7. These bright K sources at the 
center of the galaxy (L4a, L3a,b and L6 according to C05) 
are actually older clusters of more evolved stars (see Sect. 

E3D. 

NGC 5253 offers a significantly different view. It was 
already noted by Vanzi & Sauvage (|2004p and Cresci et al. 
(2005) that the radiation emitted by this galaxy is domi- 
nated at all wavelengths by one single compact source, the 
cluster identified with number 5 by Calzetti et al. (|1997|) 
at the center of our pointing. This source is a super mas- 
sive cluster with an age of just a few Myr and a mass 
of ^ 1.2 X 1O^M0, which also dominates the emission in 
our Br7 map (see Figs. [1] and H]). The Br7 source is co- 
incident with the K continuum pick. We also detect the 
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Fig. 3. Br7, [Fell] and H2 images of He 2-10, the contours of the K-band continuum are overplotted on each image. 
The coordinates are centered on the continuum peak of the brightest cluster. The crosses in the H2 images mark the 
regions discussed in Sect. 
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Fig. 4. Br7, [Fell] and H2 images of NGC 5253, the contours of the K-band continuum are over plotted on each image. 
The coordinates are centered on the continuum peak of the brightest cluster. The crosses in the H2 images mark the 
regions discussed in section 15751 



source identified as 5a by Vanzi & Sauvage (|2004p to the 
NE of source 5 and cluster 4 of Calzetti et al. p997p . Ah 
the other main clusters fall outside our field of view. Yet 
a number of fainter sources are detected, mostly in our 
continuum image without a counterpart in the ionized gas 
map, and an extended diffuse component in the emission 
line. 



3.2. Ionized iron 

The forbidden lines of [Fell] are emitted in partially ion- 
ized regions in thin shells at the edge of HII regions, or 



in more extended regions produced by shocks (Mouri et 
al. I2000p . Therefore these lines have been used as tracers 
of young supernova remnants (SNR) in starburst galax- 
ies and AGN (Moorwood & Oliva fTOMl Oliva et al 119901 
CoHna \T^M Vanzi & Rieke UMTi Alonso-Herrero et al. 
I2003p . However, the observation of these systems with high 
enough angular resolution, as well as the observation of 
star-forming regions in the Galaxy, allowed to establish 
that it is not rare to detect purely photo-excited lines in 
objects still dominated by young massive stars (e.g. Vanzi 
et al. 1^0051) . The ratio [FeH](1.64^m)/Br7 can be used 
to distinguish thermal excitation in SN shocks from pho- 
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toionization: in the first case the ratio is typicahy higher 
than a few tens, while in the second one it is significantly 
lower than 1, with typical values of 0.1 or less (Alonso- 
Herreroet aL [TM7| . 

In He 2-10 the bulk of the [Fell] emission is almost 
coincident with the Br7 one, mainly corresponding to the 
radio knots 4 and 1,2 of KJ99. However, the [Fell] peak 
is located 0.4" to the southwest of the Br7 one, and the 
emission significantly extends toward the west, reaching 
the radio source knot 3 of KJ99. A fainter source of [Fell] 
is also visible beyond the north edge of the radio source 4. 
At the position of knots 1-2 the iron is faint, while at the 
location of knot 5 only a diffuse component can be seen. 
The region with much fainter Br7 emission around knots 
1, 3, and 4 (see Sect. 13. ip is evident also in [Fell]. 

In NGC 5253 the [Fell] line is brightest at the location 
of cluster 5. It is interesting that we detect a second bright 
source 1.1" to the northwest of cluster 5, and another one 
2.6" to the east, which have no obvious counterpart in the 
continuum or in Br7. 

To better quantify the nature of these [Fell] sources, we 
used the ratio [FeII]/Br7 as a diagnostic of the excitation 
mechanism. To minimize the effect of extinction and band 
shift we used the observed ratio [FeII](1.644)/Brl2(1.641), 
and rescaled the Brl2 flux to the corresponding Br7 as- 
suming an intrinsic ratio Brl2/Br7=0.189 (Hummer & 
Storey 119871 for a case B recombination with A^ = 10"* 
cm~^ and Tg = 10^ K). The maps of the obtained ratio 
are shown in Fig. [5] The highest values are in the range 
between 1 and 7, therefore consistent with a significant 
contribution by SNe in the regions with bright [Fell] emis- 
sion. The data also support the suggestion of Johnson & 
Kobulnicky (2003) that knot 3 in He 2-10 is a non-thermal 
radio source. Indeed, this region (shown in the radio 3.6 cm 
continuum contours from Johnson & Kobulnicky (|2003p 
in Fig. [S]) is coincident with one of the regions with high 
[FeII]/Br7 ratio in our data. Moreover, C05 claim that the 
more extended radio knot 4 is a complex mix of normal 
HII regions and SNR, because it shows a non-thermal sig- 
nature with a slightly negative spectral index a. Indeed 
we find in the northern region of this radio knot the other 
region with a significantly enhanced ratio, which confirms 
their analysis. At the location of the radio knot 4 and 
1,2 in He 2-10 and cluster 5 in NGC 5253 the [FeII]/Br7 
ratio is lower, with values in the range ~ 0.1 — 0.3, in- 
dicating regions dominated by photoionization by young 
massive stars. That in both galaxies the ratio FeII/Br7 
does not reach the maximum values observed in pure SN 
remnants seems to indicate that in all cases we do not 
observe pure shock-excited [Fell] lines, but possibly an 
additional contribution from the underlying stellar pop- 
ulation. Despite this caveat, assuming that all the iron 
emission with [FeII]/Br7> 1 is dominated by the emission 
by SNR, we can estimate a SN rate for the two galaxies 
with the following relation: 



where Lip^jMtot) is the observed [Fell] flux with ratio 
[FeII]/Br7 > 1 and a nignal/noise> 3, L^p^i^{SNR) 
is the typical [Fell] emission from a SNR and T{yr) 
is the mean life of the [Fell] emission from the SNR. 
Moorwood & Oliva (fT988| . Vanzi et al. ([T997|) . Colina 
(|1993p and Mouri et al. (|2000p have computed the propor- 
tionality constant A using observations of nearby galac- 
tic SNR, finding consistent results with an average A = 
0.24 (10"*" erg/s)-^yr. Alonso-Herrero et al. (^005)) have 
instead computed the constant directly measuring the SN 
rate (SNr) in M82 and NGC 253 from radio Supernovae 
counts and the total [Fell] luminosity from HST narrow- 
band imaging. They find a lower A — 0.08, which may be 
due to a contribution from photo-excited regions to the 
total [Fell] luminosity. This value is close to the A — 0.05 
derived by van der Werf et al. (|1993p . estimating a SN 
rate from the radio synchrotron emission from the nuclei 
of NGC 6240 and measuring the [Fell] flux from the same 
regions. In this case the lower value is probably due to the 
different method to derive the SN rate. We therefore as- 
sume A = 0.24, keeping in mind that the uncertainty may 
be at least a factor of three. Using eq. ([1]) and measuring 
the [Fell] luminosity for the two galaxies in the pixels with 
S/N > 3 and [FcII]/Br7 > 1, wc find: 



SNr 



{He 2^10) 



= Ax 3.31 W^^erg/s = 7.9 W--^SN/yr{2) 



SNriSN/yr) 



L[FeII] {tot) 



L[p,ii]{SNR)T{yr) 



^xi[Fe7/](toi)(l) 



SNr(^NGC 5253) = Ax 9.68 10^' erg /s = 2.3 IQ-'' SN/yr{3) 

We compared these values with the expected number of 
SNe as derived from the far-infrared (FIR) luminosity ob- 
served in these galaxies, according to the relation between 
FIR and SN rate calibrated by Mannucci et al. ((2003)) . 
With this method, we obtain a SN rate of 1.2 x 10~^ 
and 1.1 X 10-3 SN/yr for He 2-10 and NGC 5253 respec- 
tively. The SN rates expected from the FIR luminosity 
are higher, especially for NGC 5253. However, we notice 
that our observations are not covering the whole extent of 
the activly star-forming regions of the galaxies, thus un- 
derestimating the total [Fell] flux. The result suggests a 
possible SNR origin for the observed [Fell] line emission 
in these galaxies. 

3.3. Molecular hydrogen 

The H2 line at 2.212 /ini rest wavelength and additional 
fainter lines in the H and K bands allow us to trace the 
location and physical status of the warm molecular gas in 
the two galaxies. The molecular hydrogen lines in the NIR 
are actually mainly excited through fluorescence, i.e. ab- 
sorption of UV photons, or thermally through collisions. 
The two different mechanisms mostly excite different ro- 
vibrational levels, so that it is possible to distinguish be- 
tween the two processes by comparing the ratios of some of 
the lines with models (e.g. Black & van Dishoeck ll987p . 
The 2-1/1-0 S(l) ratio is however disgnostic in the low- 
density limit, as it always thermalizes at high density (see 
Sternberg & Dalgarno [1989j) . 
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Fig. 5. Map of the [FeII]/Br7 ratio for He 2-10 (left panel) and for NGC 5253 (right panel). The radio map for He 
2-10 obtained with VLA 3.6 cm continuum image by Johnson & Kobulnicky (2003 ) is overplotted. The radio knot 3 
classified as a possible non-thermal radio source by Johnson & Kobulnicky (|2003p and the radio knot 4 classified by 
Cabanac et al. (|2005p as a mixture of normal HH and SNR corresponds with the two regions of enhanced ratio, thus 
confirming the presence of SNR driven shocks. 



He 2-10 


Line 


A (/im) 


East 


Center 


Edge 


West 


thermal 


fluoresc. 


2-1 S(3) 
2-1 S(2) 
2-1 S(l) 


2.073 
2.154 
2.248 


0.22 
0.08 
0.31 


0.21 
0.09 
0.18 


0.18 
0.07 
0.35 


0.28 
0.12 
0.27 


0.08 
0.04 
0.08 


0.35 
0.28 
0.56 


NGC5253 


Line 


A (/im) 


Cluster 5 


Cluster 5a 


SW 


W 


thermal 


fluoresc. 


2-1 S(3) 
2-1 S(2) 
2-1 S(l) 


2.073 
2.154 
2.248 


0.27 
0.22 
0.37 


0.34 

<0.13 

0.36 


0.27 
0.22 
0.43 


0.25 

< 0.09 

0.34 


0.08 
0.04 
0.08 


0.35 
0.28 
0.56 



Table 1. Comparison of the observed ratios of H2 lines to H2(l-0) SI at 2.12 /im in four different regions of both 
galaxies (see text) with fluorescent and thermal excitation models by Black & van Dishoeck (|1987p . The uncertainties 
in the ratio are on the order of 30% in He 2-10 and 15% in NGC 5253. 



As shown in Fig. [31 in He 2-10 we detect three main 
regions of H2 emission at 2.212 /im. All of them are ex- 
tended, and appear to be diffuse with filamentary compo- 
nents. They are generally not clearly associated to any fea- 
ture observed in the continuum or in other lines, except for 
the cloud corresponding to the radio knots 1 and 2, where 
Br7 and continuum emission is also detected at corre- 
sponding location. In particular, the K band nucleus does 
not show a clear counterpart in H2, and the H2 filaments 
extends beyond the bright Br7 regions, similar to what 
was observed by Vanzi et al. (|2008p in II Zw 40. However, 
if we compare our H2 observations with the study of the 
molecular gas in this galaxy as traced by the C0(J=2- 
1) transition observed by Santangelo et al. (j2009p in the 
millimeter at the Submillimeter Array (SMA) telescope, 
we find a strong spatial correlation between the CO and 



H2 clouds. Although their data cover a larger field of view 
with lower angular resolution (1.9" x 1.3"), the correlation 
between the positions of the clouds in the different species 
is clearly visible in Fig.|6l with all the detected CO clouds 
showing a counterpart in the H2. The correspondence in 
spatial distribution is not completely systematic, because 
the E-W structure to the SE of LI and the peak just east 
of L4 do not have an obvious counterpart in the CO map. 
However, a good agreement is found also in velocity, as 
shown in Sect. |3J] and Fig. [10] (see also Henry et al. l2007p . 
which provides further support for the interpretation that 
the CO and H2 we detect come from the same molecu- 
lar structures. The particular shape of the H2 emission is 
probably because we preferentially see H2 on the surface 
of the clouds, where it can be heated by UV photons or 
hit by shock waves. 
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Fig. 6. The H2 (2.212 /zm) line emission map for He 2-10 
compared with the C0(J=2-1) emission from Santangelo 
et al. (|2009p shown in contours. 



The situation is different in NGC 5253 (see Fig. g]), 
where the brightest peaks of the H2 emission can be asso- 
ciated to clusters 5, 4, and 5a. But some differences in the 
spatial extend of the emission in the ionized and molecu- 
lar hydrogen are present, because the region to the north- 
west of cluster 5 which emits in Br7 has no counterpart of 
bright H2 emission, and cluster 5 itself is elongated toward 
the southwest in H2 in a region with little Br7 and higher 
dust extinction (see Sect. 13. 4p . 

To better understand the different excitation mecha- 
nisms responsible for the H2 emission, we extracted ID 
K-band spectra from different regions in the two galaxies, 
and compared the ratios of the detected H2 lines to the 
1-0 S(l) line at 2.12 fan. For He 2-10 we extracted one 
spectrum at the location of both the eastern and western 
H2 clouds (with an aperture of 1.1" x 0.5"), as well as one 
spectrum at the emission peak at the center (0.5" x 0.5") 
and one in the gas bridge corresponding to the radio knot 
3 (1.1" X 0.25"). For NGC 5253 we extracted one spectrum 
at the peak in cluster 5 (0.5" x 0.5"), one in the cluster 5a 
to the NE (0.37" x 0.5"), and two more in the extended H2 
emission to the SW (0.37" x 0.5") and W (0.87" x 0.25") 
from cluster 5. The center of all regions is marked with a 
cross in Figs. [3] and SI In Table[T]we compare the measured 
ratios with the predictions of the Black & van Dishoeck 
(jl987p model for different excitation mechanisms. 

The 2- IS lines in the K band spectrum are particularly 
useful to distinguish between fluorescence and thermal ex- 
citation, because they are expected to be very weak in the 
latter case. Most of the lines were detected in different 
regions of both galaxies, which strongly suggests a domi- 
nant fluorescence excitation process. The measured ratios 
are lower than what is expected for pure fluorescence, al- 
lowing for some thermal contribution and suggesting that 
dense photodissociation regions are the dominant source 
of H2 emission (see also Hanson et al. I2002[l . Although the 
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Fig. 7. Extinction map in He 2-10 and NGC 5253 
from Brl2/Br7 line ratio. The K band continuum from 
SINFONI is overplotted in contours for reference. 



S/N ratio for these faint lines is not optimal for a full 2D 
spatial analysis, we obtained maps of the most significant 
line ratios and verified that there is no relevant deviation 
from the behavior sampled by the ratios reported in Table 
[TJ In particular, we do not detect any region where shocks 
alone may account for the excitation of H2, suggesting that 
dense photodissociation regions are the dominant source 
of H2 emission. 

It is an interesting point that even in the regions where 
[Fell] is strong, the H2 lines ratios does not support a dom- 
inant shock origin. This can possibly be understood by 
inspecting Fig. [51 where we see that a large fraction of the 
[FeII]/Br7 map shows values that are not strongly charac- 
teristic of shock-excited [Fell] emission (i.e. ratios in the 
range 1-2). Thus it appears that even when [Fell] emis- 
sion is strong, shock-excitation mechanisms are never ab- 
solutely dominant, although contributing to the collected 
emission. The same situation is revealed by the H2 lines 
ratios. Considering the actual nature of the regions where 
we observe the emission, i.e. galaxy centers populated by 
massive young clusters and dense HII regions, it is not 
surprising that this mixed situation is prevalent: there are 
copious amounts of UV light available as well as strong 
winds from stars and supernovae sweeping the ISM. We 
should expect to see both of these components contribut- 
ing to the excitation of emitting sources. 

3.4. Dust 

We obtained an extinction map for the two galaxies using 
the decrement of the observed ratio Brl2/Br7 with respect 
to the theoretical value of 0.189 calculated by Hummer & 
Storey fiT9~87j) for case B recombination, T=10'* K and 
n=10^ cm~^, assuming the Rieke & Lebowsky (|1985p ex- 
tinction law. The resulting extinction maps are shown in 
Fig. [7] for the two galaxies, where the contours of the K 
continuum are overplotted for reference. 

In He 2-10 the extinction is ranging between Ay ~ 
1.5 — 12, with higher values at the edges of the knots, 
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corresponding to the position of cluster L4a, and at the 
location of cluster L5 and LI. COS use the data of Vanzi & 
Rieke lTWTI to derive Ay = 10.5 from Br7/BrlO, consistent 
with the values found here. This value is also consistent 
with the ABrf—^-'2 of Henry et al. (|2007p . corresponding 
to Ay = 10.6. 

In NGC 5253 the extinction is higher at the location of 
the clusters, reaching its maximum value up to Ay = 12 in 
cluster 5. This agrees with the results of Vanzi & Sauvage 
(|2004p , who find that the brightest optical clusters of this 
galaxy are all characterized by a near-infrared excess that 
is explained by the combined effect of extinction and emis- 
sion by dust. They estimate Ay ~ 7 for cluster 5, but 
Calzetti et al. p997p concluded that the cluster has to be 
embedded in a thick gas cloud with Ay > 9 using the 
Ha/H/3 ratio, in agreement with our estimate. Moreover, 
the dust reddening is less homogeneous than in He 2-10, 
and one of the dust lanes crossing the galaxy (see Cresci 
et al. I2005II is also evident in our extinction map, south 
of cluster 5 in the NW-SE direction between cluster 5 and 
cluster 4. 



3.5. Stellar features 

The strongest stellar absorption features in the available 
spectra are the CO bands at '^ 2.3 iim, whose depth in- 
creases with decreasing stellar temperature and increas- 
ing stellar luminosity, which is therefore most evident in 
regions with red giant and supergiant stars. Following 
Leitherer et al. (|1999p . we use the spectroscopic CO in- 
dex COsp (Doyon et al. [TOMl Puxley et al. ITW|) and the 
equivalent width (EW) of the band to measure the age of 
the stellar population for the clusters old enough to show 
absorption features in their spectra. 

In He 2-10 the most prominent sources in absorp- 
tion are the cluster L4 and the arc of old clusters de- 
tected only in continuum L3/6 (see Fig. [3]). We measure a 
COsp = 0.24 and COgp = 0.22 for the cluster and the arc, 
extracting the spectra in a circular aperture of 2 pixels 
and 4x2 pixels respectively. We also measure an EW in 
the rest frame wavelength range 2.2931 < A < 2.2983 for 
the two regions of 16A and 12A. Comparing these values 
with the stellar population synthesis model Starburst99 
(Leitherer et al. rn399p , and using the Br7 detection to re- 



move the degeneracy between young and old ages for a 
given COsp or EW, we derive an age of ^ 8.3 Myr for 
the cluster LI and an age between 14.5 — 40 Myr for the 
clusters in the arc. These older clusters are lying at the 
center of the cavity between knot 1,3 and 4 detected in all 
the species discussed before (see Fig. [3] and sections 13. 1[ 
13.21 and 13. 3p . We will discuss the properties of the cavity 
and its possible origin in greter detail in Sect. 2] 

In NGC 5253 we detected CO bands at low S/N only 
for a small cluster north of cluster 5, pointing toward a 
very young stellar population in this galaxy. 
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Fig. 8. Gas kinematics from Br7 line emission in He 2-10. 
Left panel: Br7 emitting gas velocity. The velocity at the 
location of the brightest cluster is set to 0. Right panel: a 
velocity dispersion, corrected for instrumental effects. The 
contours of the Br7 line emission are shown for compari- 
son. 



3.6. Gas kinematics 

We derived radial velocities and velocity dispersions by 
fitting a function to the continuum-subtracted spectral 
profile of the emission lines at each spatial position in 
the datacube for both galaxies. The function fitted was a 
convolution of a Gaussian with one spectrally unresolved 
emission line profile of a suitable sky line, i.e. a high sig- 
nal and non-blended line in the appropriate band. A min- 
imization was performed in which the parameters of the 
Gaussian were adjusted until the convolved profile best 
matched the data. During the minimization, pixels in the 
data that consistently deviated more than 3cr from the 
average were rejected, and were not used in the analysis. 
The obtained Br7 velocity and velocity dispersion maps, 
corrected for instrumental broadening, are shown in Figs. 
[5] and O for the two galaxies. The kinematics derived with 
the same method, but using the other emission lines dis- 
cussed, are consistent with Br7, but with lower S/N. The 
errors on the radial velocities range from 5 to 20 km/s, 
and areas with larger errors are not shown in the maps. 

The gas kinematics is complex in both galaxies, with 
measured dispersions as high as 30 — 80 km/s. The high 
measured line widths are probably produced by the com- 
bination of star formation driven turbulence and outflows 
in the medium. A coherent velocity pattern is also ob- 
served in the gas kinematics of the galaxies. In particular, 
a receding area is detected in the northwestern part of He 
2-10, with a velocity difference of ^ 50 km/s with respect 
to the continuum peak. This is fully consistent with the 
CO data of Santangelo et al ((20091 see Sect. |331). In Fig. 
[TUlwe compare the C0(J=2-1) velocity map (left) and the 
H2 one (right) on the same LSR velocity scale: the two gas 
components also have a similar velocity, with the western 
cloud showing a difference of ^ 35 km/s with respect to 
the rest of the galaxy. 
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Fig. 9. Gas kinematics from Br7 line emission in 
NGC 5253. Left panel: Br7 emitting gas velocity. The ve- 
locity at the location of the brightest cluster is set to 0. 
Right panel: a velocity dispersion, corrected for instru- 
mental effects. The contours of the Br7 line emission are 
shown for comparison. 
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Fig. 10. Comparison between the C0(J=2-1) and H2 ve- 
locity field. Right panel: the C0(J=2-1) velocity map ob- 
tained by Santangelo et al. (j2009p is shown, with the 
H2 contours superimposed for reference. Left panel: the 
H2 (2.212 ^.m) velocity map from SINFONI with the 
C0(J=2-1) contours. The LSR velocity of the clouds de- 
tected in the two species agree well, with the eastern cloud 
showing a difference of ~ 35 km/s with respect to the rest 
of the galaxy. 



The molecular gas kinematics both for CO and H2 is 
fully consistent with the one traced by the ionized gas. 
This supports a picture in which the HII regions move 
at the velocities of the molecular clouds from which they 
were born. A different result was obtained with similar 
observations by Vanzi et al. (j2008p in the galaxy II Zw 
40, where at least one giant molecular cloud is detected 
with different morphology and dynamics with respect to 
the HII regions. 

We can derive a value for the V/a ratio for this galaxy 
using the observed velocity gradient and integrated disper- 
sion, finding V/a ~ 1. This value is lower than that ob- 
served in local spiral galaxies {V/a > 10), but comparable 
with the low V/a ratios observed in high-z turbulent disks 
(see e.g. Cresci et al. 12009 ', Epinat et al. 120091 Law et al. 
I2009p . Because the Jeans length ij in a gas rich system is 
proportional to {a/V)^ (Genzel et al. l2008p . this supports 
the idea that the clumpy morphology of both high redshift 
disks and blue dwarf galaxies comes from gravitational in- 
stabilities in gas with high turbulent speed compared to 
the rotation speed, and a high mass fraction compared to 
the stars (see Elmegreen et al. l2009p . 

In NGC 5253 a higher velocity dispersion area (a ^ 
50 km/s) is centered around cluster 5 (see also Vanzi et 
al. I2006p . probably produced by local episodes of winds 
or by turbulence induced by the stronger star-formation 
activity in the young central cluster. There are hints of 
an approaching area to the east of cluster 5, but the S/N 
is too low in that area to place firm constraints. No clear 
evidence of global rotation is detected in either galaxy in 
the area covered by our field of view. 

4. Feedback-induced star formation in He 2-10? 

We have already shown in Sect. l3.1ll3.2[ and l3.3l that a cav- 
ity depleted of gas is present in He 2-10 in all the species 
discussed, surrounded by some of the most actively star- 



forming regions of the galaxy, as knots 1, 3, and 4 (see 
Fig. 131). Interestingly, an arc of clusters with older stellar 
population is laying roughly at the center of this struc- 
ture, as shown in section l375l and in Fig. [3] The [FeII]/Br7 
ratio, that is a tracer of shock excited gas, is measured 
to be higher at the edges of the cavity (see Sect. 13. 2p . All 
these findings suggest a link between the formation of the 
cavity in the gas distribution and the triggering of the 
star formation activity in the surrounding regions, which 
is probably connected with the shocks responsible for the 
[Fell] emission. 

We can test this hypothesis by comparing the age of 
the clusters in the surrounding area with the expected 
travel time of a Shockwave that may be responsible for 
the formation of the cavity in the gas. Given the large un- 
certainties on both the ages of the stellar populations and 
on the velocity of the shock, a more accurate comparison 
is not possible with the available data. Assuming that the 
center of the cavity is the oldest cluster in the area L3b 
(agC'^ 12.8 Myr from Ha Equivalent Width measured by 
C05), the main clusters in the region surrounding the cav- 
ity are L4a (distance = 0.7", age-- 6.4 Myr), L4b (1.0", 
5.8 Myr), and LI (1.21", 5.2 Myr). It can be already seen 
how there is a tendency for clusters closer to the center 
to be older, as expected if an expanding shell is the trig- 
ger for their formation. If we assume a shell velocity of 
- 50 km/s (e.g. Martin [T9981 Fujita et al. [2003|) . and a 
minimum time for SNe explosions in the central cluster of 
■~ 6.3 Myr (Leitherer et al. I1999P in order to create the 
shell through feedback, we estimate that the shock shell 
reached the location of the three clusters 5.9 Myr, 5.7 Myr 
and 5.5 Myr ago, respectively. These estimates agree with 
the measured ages of the clusters, given the large uncer- 
tainties, showing that the same event might be responsible 
for the formation of the cavity sweeping away the gas and 
triggering the star formation at its edge compressing the 
denser material. 
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5. Conclusions 

We obtained integral field spectroscopy in the near- 
infrared of two nearby starburst galaxies: He 2-10 and 
NGC 5253. The study of the spectral features observed 
provided a detailed view of the environment where the 
star formation is occurring. We found that 

1. the ionized gas traced by Br7 is dominated by a single 
young super massive cluster in NGC 5253, while in 
He 2-10 the emission is due to several regions of active 
star- formation; 

2. although the ISM is mostly photo-excited, as derived 
by the [FeII]/Br7 and H2 hne ratios, some regions 
show a [FeII]/Br7 excess that is probably due to a 
significant contribution by SN driven shocks. The lo- 
cation of the [Fell] bright regions corresponds in He 2- 
10 to the positions of radio sources classified as non- 
thermal by Johnson & Kobulnicky (|2003p and Cabanac 
et al. (|2005p . We estimate the SN rates from the de- 
tected [Fell] emission, finding results that are compat- 
ible with the star-formation rates as derived by Br7 
emission and supporting this scenario; 

3. the dust extinction in NGC5253 is higher at the lo- 
cation of the main young cluster 5, reaching values as 
high as Ay = 12, shich shows that these young clusters 
are embedded in thick dust clouds; 

4. the molecular gas clouds as traced by C0(2-l) and H2 
infrared lines in He 2-10 show consistent velocity and 
morphology, and an association with the formation of 
the youngest super star cluster as traced by the ionized 
gas; 

5. we measured high turbulence in the ISM of both galax- 
ies, a ~ 30—80 km/s, driven by the high star-formation 
activity. The clumpy morphology is related to the high 
gas fraction and high gas turbulence compared to the 
rotation speed, as in high-s disks. 

6. an arc of old clusters (age > 15 — 40 Myr) is observed 
at the center of a cavity depleted of gas in He 2-10, 
surrounded by some of the most actively and young 
star-forming regions of the galaxy. This suggests a link 
between the feedback of the clusters in the arc, the for- 
mation of the gas cavity, and the triggering of the star 
formation at its edges. The estimates of the shock shell 
velocity agree with the measured ages of the young 
clusters (--5-12 Myr); 

The obtained results confirm that the unique capabil- 
ities of integral field spectroscopy can provide unprece- 
dented detail of the interstellar medium and star forma- 
tion in nearby starburst galaxies, allowing us to study the 
different mechanisms that regulates and trigger the star 
formation in these extreme environments. 
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